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Abstract 

For the case of re-jet production at next-to-next-to-leading order in the QCD cou- 
pling, in the infrared divergent corners of phase space where particles are collinear 
or soft, one must evaluate (re + l)-parton final-state one-loop amplitudes through 
0(e 2 ), where e is the dimensional regularization parameter. For the case of gluons, 
we present to all orders in e the required universal functions which describe the 
behavior of one-loop amplitudes in the soft and collinear regions of phase space. An 
explicit example is discussed for three-parton production in multi-Regge kinematics 
that has applications to the next-to-leading logarithmic corrections to the BFKL 
equation. 



*On leave of absence from I.N.F.N., Sezione di Torino, Italy. 



The quest to obtain ever increasing precision in perturbative QCD requires calcula- 
tions to higher orders in the QCD coupling, a s . Over the years, significant effort has been 
expended on computing the next-to-leading order (NLO) contributions to multi-jet rates 
within perturbative QCD. An important next step in the endeavor to obtain higher pre- 
cision would be the computation of next-to-next-to-leading-order (NNLO) contributions 
to multi-jet rates. As an example, although the NLO contributions to e + e~ — > 3 jets have 
been computed for some time now 0, the NNLO contributions have yet to be obtained. 
A calculation of these NNLO contributions would be needed to further reduce the theo- 
retical uncertainty in the determination of a s from event shape variables ||. Considerable 
effort has also been expended in the computation of leading and next-to-leading logarith- 
mic contributions to the BFKL equations for parton evolution at small x p, El. These 
two problems are, of course, connected since the logarithms that are resummed in the 
BFKL equations also appear at fixed orders of perturbation theory. Indeed, the one-loop 
Lipatov vertex may be extracted [§, 0] from one- loop five-gluon helicity amplitudes |§ . 

In order to compute a cross section at NNLO, three series of amplitudes are required 
in the squared matrix elements: a) tree-level, one-loop, and two-loop amplitudes for the 
production of n particles; b) tree-level and one-loop amplitudes for the production of 
n + l particles; c) tree-level amplitudes for the production of n + 2 particles. For the 
case of NNLO e + e~ — > 3 jets the five-parton final-state tree amplitudes, as well as 
four-parton final-state one-loop amplitudes exist in both helicity |TD[ and squared matrix- 



element forms [|TTJ, but as we discuss below, in order to be used in NNLO computations 
additional terms enter because of infrared issues. For the required two-loop three-parton 
final-state amplitudes no computations exist, as yet. Indeed, no two-loop amplitude 
computations exist for cases containing more than a single kinematic variable, except in 
the special cases of maximal supersymmetry |12| . 

Besides these amplitudes it is also important to have a detailed understanding of 
the infrared singularities that arise from virtual loops and from unresolved real emission 
when the momenta of particles become either soft or collinear. (For hadronic initial states 
infrared divergences are also associated with initial-state parton distribution functions.) 
These infrared divergences show themselves as poles in the dimensional regularization 
parameter e = (4 — D)/2. By the Kinoshita-Lee-Nauenberg theorem |13| the infrared 
singularities must cancel for sufficiently inclusive physical quantities. However, it is only 
when the various contributions are combined that the infrared singularities cancel. 

At NLO the structure of the infrared singularities has been extensively studied. The 
singularities occur in a universal way, i.e. independent of the particular particle production 
amplitude considered. Accordingly, soft singularities have been accounted for by universal 
soft functions [|TJ], [TJJ, and collinear singularities by universal splitting functions |16|]. A 
detailed discussion of the infrared singularities at NLO for e + e~ — > jets may be found, for 
example, in ref. |T7| . 
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At NNLO the situation is less developed, although some work has already been per- 
formed to illuminate the structure of infrared divergences. In particular, in the squared 
tree-level amplitudes, any two of the n + 2 produced particles can be unresolved; ac- 
cordingly the ensuing soft singularities, collinear singularities, and mixed collinear/soft 
singularities, have been accounted for by double-soft functions [I5j , double-splitting func- 
tions and soft-splitting functions [l9||, respectively. Furthermore, the universal structure 
of the coefficients of the 1/e 4 , 1/e 3 and 1/e 2 poles has also been determined [20] for the 
two-loop virtual contributions for n-particle productions. 

In this letter we shall discuss the (n + l)-parton final-state case b) when the soft or 
collinear particles are gluons. In the interference term between a one-loop amplitude for 
the production of n + 1 particles and its tree-level counterpart any one of the produced 
particles can be unresolved in the final state; the phase-space integration gives at most an 
additional double pole in e. Therefore the expansion in e of the interference term starts 
with a 1/e 4 pole, from mixed virtual/real infrared singularities, and in order to evaluate 
it to 0(e°), the (n + l)-parton one-loop amplitude needs to be evaluated to C(e 2 ). (A 
similar need to evaluate one-loop amplitudes to higher orders in e has been previously 
noted in NNLO deep inelastic scattering [^]J and in the NLL corrections to the BFKL 
equation f22|.) For the case of NNLO corrections to e + e~ — > 3 jets, this would be a rather 
formidable task given the already non-trivial analytic structure of the one-loop e + e~ — > 4 
partons helicity amplitudes presented in ref. []TU| through O(e ). 

A much more practical approach is to evaluate the amplitudes to higher order in e 
only in the infrared-divergent regions of phase-space. In the collinear and soft regions the 
amplitudes factorize into sums of products of n-parton final-state amplitudes multiplied 
by soft or collinear splitting functions. (The splitting functions in this letter are for 
amplitudes; the Altarelli-Parisi ones are roughly speaking the squares of these.) It is 
these soft or collinear splitting functions and the n-parton final-state one-loop amplitudes 
that must be evaluated to higher order in e. This is a much simpler task than evaluating 
the full (n + l)-parton final-state amplitudes beyond O(e ). 

Here we focus on the issue of supplementing one- loop (n + l)-parton final-state am- 
plitudes that are known to O(e ) with higher order in e pieces in the soft and collinear 
regions of phase space. The one-loop splitting functions have been given through O(e ) 



23], p4}| ) and the one-loop soft functions through 0(e) may be extracted from the known 



four- (^5] and five-parton || [26], 24 1 one-loop amplitudes. Below, we provide the one-loop 
gluon splitting and soft functions to all orders in e, leaving the calculational details and 
a complete listing of the one-loop splitting and soft functions, including fermions, to a 
forthcoming paper. 

As an example, we apply the framework outlined above to the one-loop amplitude for 
three-parton production in multi-Regge kinematics [27|, for which the produced partons 
are strongly ordered in rapidity. In NNLO and in next-to-leading-logarithmic corrections 



2 



to two-jet scattering, this amplitude appears in an interference term multiplied by its 
tree-level counterpart. Because of the rapidity ordering, phase-space integration does not 
yield any collinear singularities; however, the gluon which is intermediate in rapidity can 
become soft. Accordingly, the one- loop amplitude must be determined to O(e ) plus the 
contribution with the soft intermediate gluon evaluated to 0(e) p2| . To determine the soft 
gluon contribution we use our all orders in e determination of the soft functions together 
with previous all orders in e determinations of the four-gluon amplitudes 23, 29]. 



We first briefly review properties of n-gluon scattering amplitudes, since we use these 
below. The tree-level color decomposition is (see e.g. ref. |3(J for details and normaliza- 
tions) 

M* ree (l, 2, . . .n) = g {n ~ 2) £ Tr (T a ^T a ^ ■ ■ ■ T a ^) m£*(<r{l), <r(2), . . . , a(n)) , 

n 

(1) 

where S n /Z n is the set of all permutations, but with cyclic rotations removed. We have 
suppressed the dependence on the particle polarizations and momenta ki, but label each 
leg with the index i. The T ai are fundamental representation matrices for the Yang-Mills 
gauge group SU(N C ), normalized so that Tr(T a T b ) = 5 ab . The color decomposition of 
one-loop multi-gluon amplitudes with adjoint states circulating in the loop is [[H]] 

|n/2j+l 

ilC loop (l,2,...n)=^ Yl E Gr n; >)m^V(l),...,aW), (2) 

i=i (res n /s„-j 

where \_x\ denotes the greatest integer less than or equal to x, Gr„ ;1 (l) = N c Tr(T ai ■ ■ • T a ™), 
Gr n;i (l) = Tr(T ai • • • T^- 1 ) Tr(T a ^ • • -T a ") for j > 1, and S n .j is the subset of permuta- 
tions S n that leaves the trace structure Gi n] j invariant. It turns out that at one- loop the 
i can be expressed in terms of m^.} oop [0, so we need only discuss this case in this 



m n;j> 



letter. The amplitudes with fundamental fermions in the loop contain only the m^.} oop 
color structures and are scaled by a relative factor of 1/N C . 

The behavior of color-ordered one-loop amplitudes as the momenta of two color adja- 
cent legs becomes collinear, is |23, [24[ 



m 



1-loop g||& 



rv.l 



E {spht^V", b Xb ) n£3£C • • K " ■ ■ -)+SplitLl oop (a^,^) m^(. . . K\ . .)} 

(3) 

where A represents the helicity, ml l . l ° op and m^ ree are color-decomposed one-loop and tree 
sub-amplitudes with a fixed ordering of legs and a and b are consecutive in the ordering, 
with k a = zK and fc& = (1 — z)K. The splitting functions in eq. (|3|) have square-root 
singularities in the collinear limit. For the case of only gluons, the tree splitting functions 
splitting into a positive helicity gluon (with the convention that all particles are outgoing) 
is 

-1 



Split t + ree (a + ,fe+) = 0, Splitter, 



z(l — z) [a b] 
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where the remaining ones may be obtained by parity. The spinor inner products 0, |3^ 
are (ij) = and [ij] = (i + \j~), where \i^) are massless Weyl spinors of momentum 

ki, labeled with the sign of the helicity. They are antisymmetric, with norm | (ij) | = 
I [ij] I = where = 2ki ■ kj. 

The one-loop splitting functions are, 



l-loop / 



Split 



Split^ loop (a ± ,6 T ) 



{G f + G n ) Split*" 6 ! 
G n Split t i ee (a ± ,6 T ^ 



(5) 



Split 



l-loop / 



-Gf 



[ab] 



yjz{l-z) (aby 



The function G f arises from the 'factorizing' contributions and the function G n arises from 
the 'non-factorizing' ones described in ref. |34| and are given through O(e ) by [p3| , 24]] 



G f 



G h 



48vr 2 



2 e 

f + 21n(z) ln(l - z) 



(6) 



~z{\ - z)(-s ab ) 



TV 

~6 



+ 0(e) 



with Nf the number of quark flavors and 



c r 



1 r(i + e)r 2 (i- e ) 



(4vr; 



2-e 



r(l-2e) 



(7) 



As at tree-level, the remaining splitting functions can be obtained by parity. The explicit 
values were obtained by taking the limit of five-point amplitudes; the universality of these 
functions for an arbitrary number of legs was proven in ref. [54|. A listing of one- loop 



splitting functions through 0(e ) also involving fermions may be found in refs. pH|, . To 



O(e ), these splitting functions are independent of the regularization scheme parameter, 
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R 



1 HV or CDR scheme, 
FDH or DR scheme, 



(8) 



where CDR denotes the conventional dimensional regularization scheme, HV the 't Hooft- 
Veltman scheme, DR the dimensional reduction scheme, and FDH the 'four-dimensional 
helicity scheme. (For further discussions on scheme choices see refs. [^5. p5[| .) 

We have extended the above results for one-loop gluon splitting, as well as similar 
ones for soft functions, to all orders in e in several ways. The first way is by following the 



methods of ref. |34| and extending the discussion to include soft limits, but being careful 
to keep all contributions to higher order in e. In this method the contributions are divided 
into the classes of 'factorizing' contributions, that may be obtained directly from one-loop 



three-point Feynman diagram calculations and from 'non-factorizing' contributions, that 
are linked to the infrared-singular poles in e. An important ingredient in this construction 
is that the set of all possible loop integral functions that may enter into an amplitude 
are known functions to all orders in e. The method makes clear the universality of the 
splitting and soft functions since it does not rely on the computation of any particular 
amplitude. 

As a second independent method for obtaining the values of the splitting and soft 
functions we have computed the amplitudes gggH using the effective ggH coupling f36| 
due to a heavy fermion loop, again being careful to keep all higher order in e contributions. 
(For a discussion of the calculation valid through 0(e°) see ref. [0.) This is a convenient 
amplitude from which to extract the splitting and soft functions since it involves only four- 
point kinematics with one massive leg; the massive leg H ensures that the gggH amplitude 
has well defined limits when gluons are collinear or soft. As a third independent check 
we have also verified that the non-factorizing contributions obtained for the special case 
of iV = 4 supersymmetric amplitudes agree with the above determinations. (The N = A 
case has no factorizing contributions and does not provide a check of these.) The N = A 
four- five- and six-point amplitudes have been given to all orders in e in ref. |38[ making it 
straightforward to extract the collinear and soft limits in terms of limits of loop integral 
functions. 

From these calculations, our results for the all orders in e contributions to the functions 
d§) appearing in splitting functions are 



G f 



G n 



2c r 



(i-*-f)(^- 



(3-2e)(2-2e)(l-2e) v ~ " X, 
.2 \e 1 



Sab' 



(9) 



cr 



Sab, 



1 - Z 



Z 



r(l-e)r(l + e)+2 £ e k Li k ^— 



fc=l,3,5, 



—Z 



where the polylogarithms are defined as |39 

Lii(z) =-ln(l -z) 

U k (z) = f — Li fc _i(t) 
Jo t 



oo z n 



(fc = 2,3,...) 



E 

n=l 



rr 



(10) 



It is not difficult to verify that eq. @ agrees with eq. @ through O(e ). Although not 
obvious, the expression for G n in eq. (§) is symmetric in z <->• (1 — z)\ indeed its expansion 
to 0(e 2 ) may be written as 



G n 




z{l-z)\ 



1 



— + lnzln(l 



7T 



+ In z ln(l - z)-2e [Li 3 (z) + Li 3 (l 

2 



C(3)] 



+ e z 



I 2 
-lnzln(l - z)\n 2 [z(l - z)} - -ln 2 ^ln 2 (l 
6 3 



5 



+ — \nz ln(l 



360 



-7T 



+ 0(e 



The behavior of one-loop amplitudes in the soft limit is very similar to the above. As 
the momentum k of an external leg becomes soft the color-ordered one-loop amplitudes 
become 



m 



l-loop/ 
n;l \ 



.,a,k ,6, 



12) 



Soft tree (a, k ± , b) m™ p (.., a, b, ...) + Soft 1 " loop (a, k ± , 6) m^O-, a, 6, 
with the tree-level soft functions 

( ab ) Octree ^ n _ ~M 



Soft tree (a, fc+ 6) 



Soft tree (a,AT,6) 



(13) 



{ak)(kb} ' v ' ' ' [a fc] [fc 6] ' 

Following analogous methods as for the collinear case, we have computed the one-loop 
gluon soft function to all orders of e, with the result, 



Soft 1 " 10015 (a, k ± ,b) = -Soft tree (a, fc ± , b) c r 



-Sab) 



lie 



Sak)(-Skb) J sin(7re) ' 



(14) 



The soft function ( |14"D does not depend on Nf or 5r and through 0(e 2 ) it is 



Soft 1 " 10015 (a, k ± ,b) = -Soft tree (a, fc ± , b) c r 



Sab 



7T~ 



-S a k)(-S 



kb) 



7vr 4 

1 < 

6 360 



+ 0(e 3 



(15) 

Through O(e ) this agrees with the results that may be extracted from four- |2f| and 
five-parton |8|, [26], |24| one-loop amplitudes that are known through O(e ). 

We now apply these results for one-loop splitting (|5|) and soft ( |T4"D functions to the case 
of three-gluon production in multi-Regge kinematics. To do so, we also need the exact 
four-gluon one-loop amplitude through 0(e). In fact, this is known exactly to all orders 
of e. From the string-inspired decomposition of the (unrenormalized) one-loop four-gluon 
sub-amplitude p4|, |32], we write 



m 



l-loop 
4:1 



(16) 



where A 9 A , —A{, and A\ are the contributions from an iV = 4 supersymmetric multiplet, 
an N = 1 chiral multiplet, and a complex scalar, respectively. We can also write 

At = crm?»V*, x = g,f,s, (17) 

with m^ ee the corresponding tree-level subamplitude. The functions V* and V s depend 
on the helicity configuration and can be extracted to all orders in e from ref. by taking 
the massless limit. For configurations of type m^ 001 ^! - , 2~, 3 + , 4 + ) this yields, 



v' 

V s 



-h(s 23 )-e^I^ 



«13 



S23 

s 12 



e— ) I 2 D =^(s 23 ) + r» c (l - e)I, 

S12 



s 23 



D=8-2e 
4 



6 



while for configurations of type m\.{ oop (l , 2 + , 3 , 4 + ) we have 



V s 



+ 
+ 



with 



£23 r 
Sl3 



t , x . s 12 j i x S12S23 , x fD=6-2e 

MS12) H MS23) — (1 - e)/ 4 

Sl3 



'13 



•5l2 s 23( s 12 — s 23) ( fD 
'3 



„3 
6 13 



6(^ =6 - 2e (^23)-/? =,i - 2 " 



«12, 



(19) 



'13 



O12/2O23) + 833/2(312)) - — (Sl 2 / 2 D " 6 " 2£ (S23) + S 23 /- 

S13 V 

+ 



^12^23 

„3 
6 13 

S12S23 r (jD=Q-2e, 



2 D = 6 - 2e 0i 2 ; 



S12S23 ffD=6-2e 



'13 



^2*23 } D=6~2e + £12*23 ^ _ g) ^=8-2, 



'13 



'13 



Us) = ^ 



rD=G-2t, 



rD=4-2e/ 



7D=6-2e, 



JD=A-2e 
J 4 



fD=6-2e 
J 4 



fD=8-2e 
J 4 



-s/ e(l-2e)' 



-s/ 2e(l-2e)(3-2e) 

^2 ' 



1 



sy 2e(l-2e)(l-e) 

2 \ e 



l 



F l -e, -e; 1 - e; 1 + 



£23 
S12 



+ 



2(1 - 2e) 
1 

2(3 -2e) 



jD=4-2, o ;/) -: 1-2. 



vy 

D=4-2e/ 



(20) 



Fi -e, -e; 1 - e; 1 + 



£12 
S23 



+ 2I^-%s 23 ) + 2I^(s 12 ) 



£12*23^=6^ + 2 £?3^=6-2 £(s23) + 2 £l2 /^(^ 
s 13 s 13 s 13 



The functions I n are scalar loop integrals in the indicated dimension scaled by prefactors so 
as to make them dimensionless. The function V 9 obtained from the N = 4 multiplet f2"8| , 
pB] , |38|| has the same functional form for either helicity configuration. To all orders in e it 
is 

V 9 = -/f= 4 - 2e -e8 R V s . (21) 

Any partial amplitude of the type m^l° op may then be obtained from sums of permutations 
of the m^l oop pi. 



We next need the dispersive part of this amplitude in the high-energy limit, s ^> t. The 
leading color orderings of the sub-amplitudes of type m^{ oop are || (A~, A' + , B' + , B~), 
(A~,B' + ,B-,A' + ), (A-,A'+,B-,B' + ), and (A~, B~, B' + \ A' + ) where we take the Man- 
delstam variables to be s = sab, t = sbb>, and u = sab>- In eqs. ([0|) and (|2l|) orderings 
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of type mj ;1 loop (-, — , +, +) occur with su — > s and S23 — > t, while in eqs. ([H]) and (|2T| ) 
orderings of type mJ.J oop (- , +, — , +) occur with s 12 — > u and s 23 — > £ or s 12 — » i and 
s 23 - >■ M - However, for the second helicity configuration, the functions in eqs. ([T^) and 
(21) are symmetric under the exchange of S12 and S23. Thus we can limit the analysis to 
one ordering for each type. 

Using the usual prescription ln(t) = ln(— t) — iir (for t < 0), we have 



Re I 

In the high-energy limit, s 3> —t, 

t 

e, -e; 1 - e; H 



2\ e 



2 \ e 



cos(7re) . 



f22l 



2 Fi 



1 1 -orci - r) + o( — 



7T6 



Fx -e, -e; 1 - e; 1 + 



sm^7rej 
,2 \ e 



! CM — 



F a -e, 1; 1 - e; 1 + 



(23) 



S12 



El 



- + In 



S12 



+ 



with S12 -> s or S12 — > it. Since u = —s — t, for either choice of Si 2 the dispersive parts 
are the same to the required accuracy, and we can take Si 2 — > s in eqs. (p3j). Substituting 
eq. (|22"D and (|23|) in eq. (|20|) , and using the identity 



cos 7re 



7T- 



sm 7re 



V(l + e)-V(-e) 



(24) 



we obtain 



Disp 7, 



D=4-2e 



2 //x 



2 \ e r 



^(l + e)-2^(-e)+V(l)+ln- 



+ O 



(25) 



In addition, in the high-energy limit both eq. (|IBD with S12 — > s and S23 — > t, and eq. fll9| ) 
with S\2 — ¥ u and s 2 3 — > t reduce to 



v' 

V s 



1 

-t / e(l - 2e) 
V 1 



of* 



(26) 



e(l - 2e)(3 - 2e) 

while the function V 9 , eq. (pT|) , is obtained from eq. (p5|) and (|26|). 

Using eqs. (|16|), (0), (HI), (|25|) and (p6j), and the fact that the proportionality factor 
between each tree-level subamplitude and its one-loop correction is the same for all color 



S 



orderings, we obtain the unrenormalized four-gluon one-loop amplitude in the high energy 
limit to all orders in e, 

,2 \ e 



DispM!" 1 °° V {A~, A' + , B' + , B~) = Ml me (A~, A' + , B' + , B~ 



g 2 c r 



x 



2(1 - 2e) ^(1 + e) - 2 ^(-e) + VK 1 ) + In — ) + 



1 



el 



3-2e 



4 



+ 



2c) 
2(1 



(27) 



3-2e 



-iVy 



Following the methods of ref. ||, we can then extract the unrenormalized one-loop cor- 
rection to the helicity- conserving vertex, to NLL accuracy, 



DispCif\-p a , Pa ,) DispC£ ! 



99(1)/ 



-Pb,Pv 



C 



■99(0) . 



'Pa, Pa' 



(1 - 26) 



-y99(0) , 



c r 



1 



1 



C-+ '{-p b ,Pb> 
- e)-2^(-e) + </>(!)]- 



-t J e(l - 2e) 



(2* 



- 2 



1 



3-2e" 



/ 



2(3 - 2e) 

Eq. ( EH ) is valid to all orders in e for <5r = and 1; it agrees with the one- loop correction 
to the helicity-conserving vertex computed in ref. || to O(e ) and with the one computed 
in ref. []27|, [40| to all orders in e, for 5r — 1. 

In an inclusive high-energy two-jet cross section at NNLO and in the NLL corrections 
to the BFKL kernel, the five-gluon one-loop amplitude is multiplied by the correspond- 
ing tree-level amplitude with the intermediate gluon k integrated over its phase space. 
Assuming multi-Regge kinematics, the only infrared divergence that can arise is in the 
soft limit for the intermediate gluon, k — > 0. It gives a single pole in e. Thus, in order 
to generate correctly all the finite terms in the squared amplitude, the five-gluon one- 
loop amplitude must be computed exactly to O(e ), and must be augmented by the (9(e) 
corrections in the soft limit for the intermediate gluon. To achieve that, we need the 
dispersive part of eq. ( fL2l) in the physical region, s a b > 0, s a k > 0, Skb > 0. (The other 
leading color orderings yield the same result.) Using eq. flllD , (^2]), and the identity ([24]), 
we can write the dispersive part of the soft function to all orders in e as, 



Disp Soft 1_loop (a, fc ± , b) = -Soft tree (a, k ± , b) 



cr 



In addition, using the strong rapidity ordering, 



y a ~>y~> w> 



\k 



a±\ 



' V? S ab " 
Sak Skb , 



[l + e^(l 



I*. 



6_L 



the mass-shell condition for the intermediate gluon gives 

Sak Skb 



Sab 



Subsequently, eq. (|29f) becomes, 



Disp Soft 



l-loop . 



I A; i 



0,^,6) = -Soft trce (a, hT,b)-j 



-e^(l + e)] . 

(29) 

(30) 
(31) 

e^(l + e)] . 

(32) 
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The soft limit for the tree-level five-gluon sub-amplitudes, 



tree 



{AT, A' + , fc* B'+, B~) = Soft tree (A', k ± , B') mfCA" A'+ £'+ £" 



(33) 



holds for arbitrary kinematics. Thus, the unrenormalized five-gluon one-loop amplitude 
in the multi-Regge kinematics, and in the soft limit for the intermediate gluon and to all 
orders in e, is obtained by using eq. fll2|) , with the four-gluon one-loop amplitude 
the loop soft function (0), and eq. ([3~3~D , yielding 



Disp M 5 1-loop (A-, A' + , k±, B'+, B 



g 2 c r M 5 tree (A", A'+ fc* B' + , B~ 



To 0(e), eq. 




+ ^ ^(1 + e) - 2^(1 - e) + ^(1) + In 



(34) 



1 - S R e 



2e 



+ 



2(1 



e(l - 2e)(3 - 2e) 



- [l + e^(l-e)-e^(l + e)] 



DispM 5 1 " loop (A~, fc* B'+ S- 



c/ 2 c r M 5 tree (A", A'+ fc* B' + , B' 



2 \ € 




+ - In 



-t 



A) , Ar /10 56 

h i\f \ e 

e H 9 27 



\k\ 



-y + 2C(3)e -^r 



+ 0(e 2 



9 



5*e (35) 




with (3 Q = (11N C - 2N f )/3. We have checked that eq. (gD agrees to O(e ) with the 
five-gluon one-loop amplitude, with strong rapidity ordering and in the soft limit for the 
intermediate gluon || |7|. The above result may be used to verify the virtual next-to- 
leading log corrections to the Lipatov vertex for use in the BFKL equation as is done in 
ref. §]. 

The same type of analysis may be applied more generally to the problem of NNLO 
QCD corrections. The one-loop gluon splitting and soft functions that we have presented 
here are valid to all orders in the dimensional regularization parameter, e. This allows 
them to be used in NNLO calculations with infrared singular phase space where terms 
of up to two powers in e are necessary. Previous explicit determinations of the one- 
loop collinear splitting functions |23|, 24| were not performed to the required order in e. A 
systematic discussion of the soft and collinear splitting functions and further calculational 
details, including the case of external fermions, will be presented elsewhere. In particular, 
these functions can be used to aid in the computation of NNLO contributions to e + e _ — > 3 
jets once all the matrix elements are available. However, much more remains to be done 
before this is realized. 
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